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Abstract 

In this paper, the isothermal plasma wave properties in the neigh- 
borhood of the pair production region for the Kerr black hole magneto- 
sphere are discussed. We have considered the Fourier analyzed form of 
the perturbed general relativistic magnetohydrodynamical equations 
whose determinant leads to a dispersion relation. For the special sce- 
nario, the x-component of the complex wave vectors are numerically 
calculated. Respective components of the propagation vector, atten- 
uation vector, phase and group velocities are shown in graphs. We 
have particularly investigated the existence of a Veselago medium and 
wave behavior (modes of waves dispersion). 
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1 Introduction 

The sky is full of fascinating cosmic objects including planets, stars and 
compact objects like black holes. A black hole is a region of space where 
the pull of gravity is so strong the even light cannot escape. Black holes are 
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described in the framework of relativistic astrophysics. General Relativity 
leads to insight and understanding of black hole physics. 

Rotating black holes exist in the nuclei of galaxies. The gravitational 
pull of the black hole drags the material (highly magnetized plasma) from 
nearby stars and forms an accretion disk as it spirals inwards and finally 
disappears into the black hole. Highly magnetized plasma surrounding the 
black hole horizon is considered as a magneto-rotational fluid governed by 
the general relativistic magnetohydrodynamical (GRMHD) equations due 
to its motion under immense gravity. The magnetohydrodynamical waves 
produced in these magneto-rotational fluids not only help to confirm the 
existence of a black hole but also transmit the information inside the magne- 
tosphere. In addition, these show the response of black holes to the external 
perturbations. The 3-1-1 ADM formalism [1] is frequently used to discuss the 
space plasmas into the arena of black holes. Thorne and Macdonald [2]-[l] 
dealt the electromagnetic fields of the black hole theory with 3-1-1 formalism. 
Later, this formalism was used by Holcomb and Tajima [5], Holcomb [6] and 
Dettmann et al. [7] to study wave propagation in the early universe. Buzzi 
et al. 011] used the 3-1-1 GRMHD equations to study one dimensional trans- 
verse and longitudinal waves in two-component plasma in the vicinity of the 
Schwarzschild horizon. 

There exists large body of literature [10]-[27] which indicates keen interest 
in plasmas, plasma processes and plasma waves for the magnetospheres of 
compact objects. Punsley fnj\ discussed the electromagnetism and plasma 
waves created in the ingoing and outgoing perfect MHD ergospheric winds. 
Beskin et al. [11] gave the comprehensive analysis of the pulsar magneto- 
sphere and discussed the pair creation region around a black hole. Park and 
Vishniac [12] included the effects of plasma accretion flows on the mass and 
angular momentum-loss rates in addition to the pure electromagnetic extrac- 
tion of energy from the rotating black hole. Takahashi et al. [13] focused on 
the Alfven surface. By considering the cold limit, trans- Alfvenic solutions 
are characterized and MHD inflow is discussed in relation to the extraction 
of energy from the rotating black hole and the accretion of winds or jets. 

Ruffini and Wilson [14J proposed the idea of energy extraction from the 
Kerr black hole. Blandford and Znajek [15] were the first to construct a global 
model for the magnetosphere of a Kerr black hole in the force-free limit. They 
demonstrated that energy could be extracted under certain conditions in the 
form of a Poynting flux. They applied their model to radio loud AGNs, and 
proposed that extragalactic jets are the consequence of this energy extraction 
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mechanism [TB] . Several authors [T7j discussed different aspects of Blandford 
- Znajek process. 

Energy extraction is discussed by different authors pE] assuming the mag- 
netic field lines connected with the transition region instead of connecting 
with remote loads inside the accretion disk. Zhang [191 120] formulated the 
3+1 black hole theory for stationary symmetric GRMHD and applied it to 
cold plasma filled Kerr magnetosphere. Using the interface conditions defined 
at the pair production region, he obtained the numerical solutions which show 
that outflow of energy flux is possible for the specific angular frequency and 
specific x-component of the wave vector. 

Lee et al. [2T]-[23] considered the Blandford-Znajek process as one of the 
viable models of powering the gamma-ray bursts effectively. The gamma ray 
bursts is a consequence of interaction of magnetic field and plasma. The 
observations of a broad Fe Ka line in the bright Seyfert 1 galaxy MCG-6-30- 
15 and M87 plj-[26] suggest evidences of the extraction of rotational energy 
from a Kerr black hole by a magnetic field. Thus, the popularity of the 
GRMHD behavior of plasma has increased and the study of GRMHD waves 
become more important. Koide et al. [27] discussed plasma flowing into a 
rapidly rotating black hole by using numerical simulations via GRMHD. 

Sharif et al. [28] considered the cold and isothermal plasmas in the vicin- 
ity of a Schwarzschild event horizon. Plasma wave properties are investi- 
gated using 3+1 formalism. The results are generalized by considering the 
Schwarzschild geometry as a restricted Kerr model. We [2n]-|2I] have inves- 
tigated cold plasma wave properties (real and complex wave numbers) and 
isothermal plasma wave properties (only real wave numbers) |32] for Kerr 
planar analogue. 

This paper is devoted to discuss the isothermal plasma wave properties 
by using complex wave numbers. We consider an isothermal plasma in the 
vicinity of the Kerr magnetosphere's pair production region. The correspond- 
ing dispersion relations are obtained by using Fourier analyzed perturbed 
GRMHD equations. We then evaluate wave numbers which are used to find 
x-components of propagation and attenuation vectors, phase and group ve- 
locities. These quantities help to investigate the properties of plasma. 

The organization of the paper is given as follows: Next section provides 
the background metric with flow assumptions and respective physical inter- 
pretation. Section 3 contains the Fourier analyzed GRMHD equations for 
isothermal plasma in Kerr planar analogue. Section 4 is devoted to the 
numerical solutions of the dispersion relation including Figures and related 
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discussion. We shall summarize this discussion in Section 5. 

2 3+1 Background and Relative GRMHD Equa- 
tions 

The general line element in ADM 3+1 formalism can be written as [20] 

ds^ = -a^dt^ + riij{dx' + p'dt){dx^ + p^dt), (2.1) 

where a is the time lapse and (3 is the three-dimensional shift vector and rjij 
are the components of three-dimensional hypersurface metric. The quantities 
a, Pi, rjij, = 1,2,3) are dependent on coordinates (t, Xi). 

We consider the planar analogue of the Kerr metric (Eq.(2.2) of [32J) 

ds"^ = -dt^ + {dx + P{z)dtf + dy"^ + dz^ . (2.2) 

Here the lapse function is set to unity and the shift vector P is along x- 
direction only. The details of the analogue metric are mentioned in Section 2 
of [32] • In this planar analogue (Eq. fl2.2p ). we consider a rotating isothermal 
plasma admitting the following equation of state 

jjL = ^ ^ = constant, (2.3) 
Po 

where /x, po, p and p (all are taken to be constants) are the specific enthalpy, 
rest and moving mass densities and pressure of the fluid respectively. This 
highly magnetized plasma includes the effects of gravity (due to nearest cor- 
respondence to the black hole). Thus perfect GRMHD equations express it 
completely. 

For the line element given by Eq. fl2.1l) . the perfect GRMHD equations 
(Maxwell's equations of magnetic field evolution, local conservation laws of 
mass, momentum and energy in force-free magnetosphere respectively) are 

m 

^ + -B.V/3 + ^B = -Vx (aVxB), (2.4) 

dr a a 
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■ Po-y'^.^ + «i I + v.(«v - /?) H 0. (2.5) 
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^ -(B X {V X [V X (aV X B) - (B.V) 



+ (V X B).V^})„ (2.6) 
|:(po/i7^) - ^ + ^ ((PoP7^ - + i-((V X B)2 + B^)) 

+— {PoP7^V^'V^^' - —((V X B)''(V X B)^' + B'B^) 

2a An 

+ ^((V X B)2 + B^)7]'^}£tVi, + Po/^7'V.a + po/i7'V.(V.V)V 

OTT 

_±a.((V X B) X B) - -^(V X (V X B)).B - -{po/i7'V.(V.V)^ 

ZTT 47r a 

+p{^.P) - -^(V X B).((V X B).V)^ - -^B.(B.V)^ 

477 An 

+ ^((V X B)2 + B')){V.m + ^(-{Vx B).{(V X B).V}^ 

+ (V X B).(a X B) + 2(V x B).^(V x B)(3 

dr 

dB \ 

+^(VxB).(VxB) + B.— =0, (2.7) 

dr J 

where B, V, a respectively denote the three dimensional magnetic field, ve- 
locity and acceleration of the fluid, 9 gives the expansion of the fluid and ^ 
is time derivative measured by the flducial observer, is the time deriva- 
tive with respect to motion of the fluid, £t denotes the Lie derivative with 
respect to time parameter t. The quantity g is the determinant of the three 
dimensional hypersurface metric. The respective form of Eqs. (l2.4p - fl2.7l) for 
the Kerr geometry is given in Appendix A. 
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3 Plasma Discussion and Physical Interpre- 
tation 



We are going to perform a local analysis to the plasma waves. For this 
purpose, we have assumed that the plasma living in the spacetime fl2.2l) 
admits a free fall along z-direction due to the black hole's gravity and along 
x-direction due to rotation of the black hole (in direction of shift vector of 
our analogue spacetime). The fluid's velocity (V) and magnetic field (B) 
measured by our fiducial observer (a natural observer associated with the 
spacetime and locally measures the quantities) can be described as follows: 

Y = V{z)e^ + u{z)e,, B = 5{A(^)e, + ej, (3.1) 

where i? is a constant. 

In our stationary symmetric background, the x-component of the velocity 
vector can be written in the form [20j V = C + Am, where C = (3 + Vf with 
Vp as an integration constant. We assume the value of the shift function [20] 
13 = tanh(2;) — 1 and Vp = 1. Further, = Stt and A = 1 are taken so that 
the magnetic field becomes constant. Thus the x-component takes the form 
V = 1 + P + u = tanh(2;) + u. Substituting the value of V in the conservation 
law of rest-mass for three-dimensional hypersurface, i.e., 

Po'ju = fi{p + p)'ju = A (constant) (3.2) 

with the assumption that rest-mass density is constant and A/ po = 1, we get 
an equation of the form 

3m^ + 2Mtanh(z) + tanh^(2) -1 = 

quadratic in u. The solutions of this equation [ui and U2) and the corre- 
sponding values of V {Vi and V2) are given as follows (Eqs.(5.2)-(5.3) of 
i2): 



ui = -|tanh(z) - 1^3 - 2tanh^(2;), I 

Vi = |tanh(2) - |y^3 - 2tanh2(2), 

U2 = — I tanh(2;) + |^3 — 2 ta,nh^{z), 

V2 = |tanh(z) + 1^3 -2ta-n.h\z). 
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The corresponding Poynting vector (S) is 

S = -;^ExB = 2 tanh(^) (ex - e^) , (3.5) 

where E is the electric field as measured by the fiducial observer. This electric 
field can be determined by the perfect fiow assumption, i.e., 

E = -V X B. (3.6) 

We have assumed that plasma is produced in the pair production region 
lying at z — 0. A perturbation method is used in which we perturb the 
equilibrium by a very small amount, for example, the background plasma 
mass density is not the same as po everywhere but changes slightly with 6p. 
The perturbed quantities are assumed to be small and we have considered 
these as first order quantities. Perturbation assumptions and equations are 
given in Appendix A. 

The perturbations are produced due to plasma production in the pair 
production region. In this region, electron-positron pairs are produced and 
taken apart by strong forces of gravity. The electrons move towards the 
outer end of the magnetosphere while the positron move towards the event 
horizon. Thus the displaced electron results a positive charge density which is 
produced at the position from where electrons are displaced. These positive 
charges then attract the electrons. Electrons move back but they overshoot 
it, and come back but overshoot again and thus oscillation is produced. Same 
is the case of positrons which produce oscillations in the region towards the 
event horizon of the black hole. 

Here we check the properties of the plasma as a consequence of the plasma 
perturbations. For this purpose, we have specifically tried to over check the 
regions in the neighborhood of the pair production region, i.e. z = 0. It is 
important to note that the production of plasma does not make any change 
in the background geometry. It only disturbs the fiuid dynamics. Thus, all 
geometric terms are considered to 0th order. 

4 Fourier Analyzed Perturbed GRMHD Equa- 
tions 

We have assumed that the perturbations produced are simple harmonic 
waves. We shall discuss the wave properties with the normal mode anal- 
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ysis (Fourier decomposition of waves). The following assumptions will be 
used for the sake of Fourier analysis (Eq.(3.3) of [32]). 

Pit, X, z) = c^e"^^'-^^^-^^''\ pit, X, z) = CQe-'^^t-k.^-k.z) ^ 
v^it, X, z) = c2e-*("*-'=-"-'=^"), v,it, X, z) = cge-'^"*-'^-"-'^-"), 
hit, X, z) = c4e-*("*-^-"-^'^"), hit, X, z) = c5e-'("*-'=-"-'=^"). (4.1) 

Here q, (i = 1, 2, 6) are constants, k = (fej,, 0, k^) is the wave vector and 
u is the angular frequency of the waves. 

The wave vector leads to the following quantities: 

• i?e(k) gives a vector, called propagation vector, that shows three di- 
mensional propagation of waves. 

• Jm(k) gives a vector, called attenuation vector, that shows damping 
and growth of waves in three dimensions. 

• The phase velocity vector gives the speed of phase in three dimensions 
and can be calculated by the formula 

• The group velocity vector gives the speed of wave packet in three di- 
mensions and can be calculated by the formula ^^j^-^ . 

We consider the perturbed perfect GRMHD equations written in 3+1 
formalism for the Kerr planar analogue with isothermal state of plasma 
(Eq.f l2.3p ) given by Eqs. (lA13l) - (lA14l) . The Fourier analyzed perturbed form 
of these equations is given as follows (Eqs.(4.6)-(4.12) of [32]): 

LkzC2 — iik^X + A')c3 — c^iikzU — lu + u') + c^i'dikz + = 0, (4.2) 

LkxC2 — ikxXc^ + LCK,i'dkx + kzU — cj) = 0, (4.3) 

/c^C4 = -/C2C5, (4.4) 
Ci[Lpi-uj + §k.j, + ukz) - ip'u + pu' + pM7^v^)] + C2ip + p) 
X [-r/7V + Lk,-^^uV + Lk^-K + -i'^u{i'K + -i^V^)V' + 



+C3(P + P) 



u 

-ri-f'^u + Lk^j'^uV + Uk, - (1 - 27^^^)-)^ + 

u 

,,2, 



-\-CQ[ipi—uj + -dkr^ + ukz) + ip'u + pu' + pw) (f)] = 0, (4.5) 
Cip'y'^u{7rV' + -f'^uVu'} + Cqp[-i'^u{W + '-f'^uVu'} + ik^] 
+C2 H{(P + P)7'7r + ^} +C{(P + P)7'(1 + 7V2) -V^} 



+ (P + p)7^M{(7r + 37V2)W + AttVV'}] + Cg [-7] {{p + p)^\V 
-iPX} + C {{p + ph^uV + + {p + p)^^[2-f\Q + -f\^)uVu' 
+{{g + j\^)in + 7V) - -f^V^}V'] + i:u\'] + i:c^{-tK{l - u^) 
+uu'} + ipC5 {-A' - wd' + - 1/2) - 2LuVk,} = 0, (4.6) 
ci-f'^p[u{gu' + -f^VuV'} - V(3'] + ce[-r'^p{guu' + 7VmV' - Vf3'} 
+p' + ptk;,] + C2 [-7] { (p + p)'y^uV - ^X} + C { (P + p)l'^uV 

+ (p + p)72{7 W(7r + 37V2) - /?'(7r + 7V2) + 2V-f^uu' 
x(f? + 7V)}] +C3 + + +C{(p + P)7'f?-^A'} 

+ (P + P)7'K(1 + 7'«')(1 + 47'n') + 2M7'{(f? + j\^)VV' - VP'}] 
-ipXX'u] + ipc^ {ik^X{l - u^) + A' - AW} + i)Cr^{2XuVLk^ + Am^?' 
-AtA;,(l-y2)] =0, (4.7) 
Ci7^[p'u — p(iCo' + u' + 2u7^(y9 — 7^u\//3' + Lkx'd + ifc^w)] 
+C6[-iu;p(7^ - 1) + I'^ip'u + p(m' + 2M7V - l'^uVl3' + ifc^-i? + Lk^^u)}] 
+C2[-Luj{2{p + p)7V - ^x} + ^A;.[(P + P)7'{1 + 2^^Vd} - ^jdx] 
+ik,u{2{p + p)7V + ^/-x} + (P + pWu{2^^V' + 67VV9 - /3'(7r 
+7V2)} - ^A'] + C3[-6cj{2(p + p)7^M + ^Ax} + iA;,i9{2(p + p)7^M 

+7/'Ax} + ik^{{,P + p)7^(^ + I'^u'^) - ^Amx} + (p + P)7^{ + 27^^' 

+67V¥P + 7V - + V^A'{A - Mx}] + C4V^[m{A' - xu] + ik,x 

X (1 - u^)] + C5^{-AV + Mx^' - 6fc^x(l - ■^^^) + 2Lk,uV} = 0. (4.8) 

where 7^ = V - f3, = W + mm', x = uX-V, V' = £, ^1 = 1 + 7^^^^ vr = 
1 + 7^^^, C = /-{kxV + kzu) and rj = t{uj + pk^). Here prime denotes the 
derivative of the quantity with respect to the variable z. 



5 Numerical Solutions 

Equation (14. 4p gives k^ = ~^kx and we will assume k^ = —k^ for the sake 
of simplicity. Using this relation and the assumptions given in paragraphs 
1 and 2 of Section 3, we simplify Eqs.(i2D, (Si. (gSD-dSHD- To obtain a 
dispersion relation, we use the method given in [33]. We equate the deter- 
minant of the coefficients of constants (ci, C2, C3, C4, C5, cg) for Eqs. (l4.2p . 
(14.31) . (I4.5p - (I4.8I) to zero. It gives an equation of the type 
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+A5{uj, z)kj + Ae{uj, z)k^ + Aj{uj, z) = 0, (5.1) 



sextic in k^ which cannot be solved exactly to obtain the values of k^. Thus 
we use the software Mathematica to solve it by using the fluid flow velocities 
(13.31 and l3.4p . Related Mathematica codes are given in Appendix B. 

We assume that the pair production region is exactly at z = 0. Our region 
of consideration is —5 < z < 5 for wave analysis. This region is taken only 
to discuss the waves near the pair production region as much as possible. 
Since the flow variables have large variations in the region — 1 < ^ < 1, 
we avoid this region and solve the dispersion relation for —5 < z < —1 
and 1 < z < 5 and thus we have two meshes (just like the near and far 
zones in em-radiation). The mesh —5 < z < —1 indicates the neighborhood 
of the pair production region towards the event horizon and the mesh 1 < 
z < 5 represents the neighborhood of the pair production region towards the 
outer end of the magnetosphere. We deal with each mesh separately. For 
a numerical solution, we take the step-length 0.2 for z and u and find the 
value of kx at each point {zi,ujj), {i = 1, ...,21, j = 1, ...51) of the mesh. 

Our dispersion relation (Eq. (l5.9p ) leads to six complex values of kx at 
each point. Thus k^ = k^^ + ^^/x- The real part kRx represents the x- 
component of the propagation vector which leads to the x-components of 
the phase and group velocities. The imaginary part kj^ represents the x- 
component of the attenuation vector. The x-component of the propagation 
vector gives the x-components of the phase velocity {vpx = -j^) and group 
velocity {vgx = -^f^) of the waves. Each root of the dispersion relation is 
separated and numerical interpolation is used to estimate the surface. The x- 
components of the propagation vector, attenuation vector, phase and group 
velocities can be evaluated from these interpolation functions. 

It is mentioned here that we investigate the x-component of the above 
quantities. The equation k^ = —k^ indicates that the behavior of the z- 
component of the wave vector is opposite to that of the x-component. If 
the x-component of the wave is normally dispersed in one region, the z- 
component is anomalously dispersed and vice versa. When the x-components 
of the propagation vector, attenuation vector, phase and group velocity vec- 
tors increase, their corresponding z-components show a decrease and vice 
versa. The components of the propagation and attenuation vectors, the 
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phase and group velocities lead us to the wave behavior of the Kerr black 
hole magnetosphere and the properties of a Veselago medium. 

Here, we give a criteria which is used to discuss the two properties of the 
medium. 

• A medium is a Veselago if the propagation vector is in opposite direction 
to the Poynting vector [M]. In the usual medium, both these vectors 
admit same direction. 

• If the phase velocity is greater than the group velocity, dispersion of 
waves is normal, anomalous otherwise |35j . 

We shall use these facts to infer our results. 

We have considered two velocities (given by Eqs. (13.31) and (13.41) ) for the 
rotating plasma. For each velocity, we obtain six complex values of the x- 
component of the wave vector. For fluid flow with velocity components given 
by (13. 3p . Figures 1-6 show the dispersion relations in the neighborhood of 
the pair production region towards the horizon whereas Figures 7-12 indicate 
these in the neighborhood of the pair production region towards the outer 
end of the magnetosphere. For the velocity components given by Eq. (l3.4p . 
Figures 13-18 and 19-24, respectively show the dispersion relations in the 
neighborhood of the pair production region towards the event horizon and to- 
wards the outer end of the magnetosphere. All the Figures have graphs A, B, 
C and D representing the x-components of the propagation and attenuation 
vectors, phase and group velocities respectively. 

In the following, we shall discuss some of the figures which lead to im- 
portant results. A summary of the discussion of results obtained from the 
remaining figures is given in the last section. 

In Figure lA, the x-component of the propagation vector takes nega- 
tive values in the regions —5 < z < —4.598 and —4.4 < z < —1 while 
it has positive values for —4.598 < z < —4.4. For the region towards the 
event horizon, the x-component of the Poynting vector, given by Eq. (l3.5p . 
admits negative values. When we compare these two, the x-components of 
the propagation and Poynting vectors turn out to be in the opposite direc- 
tion for the region —4.598 < z < —4.4 and thus the medium is a Vese- 
lago. Figure IB indicates that the x-component of the attenuation vector 
increases with the increase in angular frequency and decreases in z in the 
region —2.0 < z < —1.0, 0.175 < u; < 10. This shows damping of waves 
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with increase in angular frequency and wave growth with increase in z. The 
remaining region shows random damping and growth of waves. 

Similarly, in Figure 7A, the x-component of the propagation vector is 
negative throughout the region. The x-component of the Poynting vector, 
given by Eq. fl3.5p . admits positive values for the region away from the event 
horizon. The opposite direction of the x-components of the propagation and 
the Poynting vectors indicates the presence of a Veselago medium. 

Table 1 gives the figures along with the region where we have found a 
Veselago medium, i.e., the the propagation vector is in opposite direction to 
the attenuation vector. 

The negative phase velocity propagation is a property of the usual medium 
as well as Veselago medium. In this medium, the right hand rule of electro- 
magnetism changes into the left hand rule. This medium has been demon- 
strated experimentally |36] (the material formed is called left-handed material 
or meta-material) and is characterized by simultaneous negative magnetic 
permeability and electric permittivity. It has Poynting vector of a time- 
harmonic plane wave directed opposite with respect to the wave vector. This 
exhibits unusual properties like support of backward waves and anomalous 
negative refraction of plane monochromatic electromagnetic waves. Thus, 
backward waves are supported by the plasma in these regions. 

Figures IC and ID indicate that the x-component of the group velocity 
is greater than the phase velocity in the region —1.75 < z < —1.0, 0.5 < 
C(j < 10 which shows anomalous dispersion of waves. Rest of the region 
admits random points of normal and anomalous dispersion of waves. Table 2 
contains figures with regions of the normal dispersion of waves. The regions of 
normal dispersion indicates that the waves can pass through them without 
any restriction. Thus the energy flux can pass through these regions (via 
waves) . 

6 Summary 

This paper is devoted to discuss the properties of isothermal plasma waves by 
using the perturbed GRMHD equations in 3+1 formalism. These non- linear 
ordinary differential equations are simplified by using the Fourier decompo- 
sition of waves (Fourier analysis of waves). We have computed the dispersion 
relation for the sinusoidal waves and solved it to check the properties of the 
medium. This has been done in the neighborhood of the pair production 
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Figure 


Regions of Veselago Medium 


2A 


-1.8 <z<-l 


3A 


-2.05 <z< -1.0 


4A 


-2.1 <z< -1.0 


5A 


-5<z<-l 


6A 


-5<z< -4.6 and -4.375 < z < -1 


8A 


l<z<5 


14A, 15A, 16A, 17A and 18A 


-5< z< -1 


19A 


0<z< 4.425 and 4.585 < ^ < 5 


24A 


4.42 <z< 4.595 



Table 1: Regions admitting a Veselago medium in respective figures 



Figure 


Regions of Normal Dispersion 


2C and 2D 


-1.8 <z< -1.0 


7C and 7D 


l<z<5 


8C and 8D 


0.5 < < 5 


IOC and lOD 


0.15 <uj < 10 


12C and 12D 


l<z<5, 0.3 < a; < 0.401 


13C and 13D 


0.1 < a; < 0.5 


14C and 14D 


< a; < 0.5 


15C and 15D 


Except for negUgiblc angular frequency 


16C and 16D 


-5.0 <z< -1.025, 0.4 < a; < 10 


18C and 18D 


1.01 <uj <2 


19C and 19D 


l<z<2, 0<a;<10and2<z< 2.8, 0.75 < a; < 10 


20C and 20D 


l<z< 1.8, 0.5 < a; < 10 



Table 2: Regions of normal dispersion of waves in respective figures 
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region of the Kerr magnetosphere. 

A dispersion relation is found by using MathemafAca and wave numbers 
are deduced in terms of angular frequency. The real part of the x-component 
of the wave vector gives the x-component of the phase and group veloci- 
ties which lead to the properties of plasma. The imaginary part of the x- 
component of the wave vector indicates the x-component of the attenuation 
vector. All these quantities are shown in the graphs. 

The following two objectives are investigated in this work: 

1. The behavior of GRMHD waves under the influence of gravity and 
magnetospheric wind is analyzed (especially dispersion of waves is con- 
sidered). This helps us to detect whether the extraction of energy is 
possible from the black hole. 

2. The existence of a Veselago medium in the black hole regime is checked. 

To obtain these objectives, we considered two values of the flow veloc- 
ity. The dispersion relations for each velocity are calculated and complex x- 
component of the wave vectors are evaluated by using the relation k,j. = —k^. 
This relation gives us the ;2-component of the wave vector dependent upon 
the x-component. The direction of the Poynting and the wave vectors in 
the region gives an estimate of the refractive index (positive if both have the 
same direction and negative if both have the opposite direction). The effect 
of rotation of the black hole yields and the gravitational effect gives kz- 
The results can be summarized as follows: 

In Figure 1, the medium is usual near the pair production region and 
shows anomalous dispersion for most of the part. Figure 2 indicates that 
normal dispersion of waves occurs in a Veselago medium near the pair pro- 
duction region. Figure 3 shows random points of normal and anomalous 
dispersion in a Veselago medium living near the pair production region. Fig- 
ure 4 shows that there exists a Veselago medium with random normal and 
anomalous points of dispersion. In Figure 5, a Veselago medium is found 
with random points of normal and anomalous dispersion in the neighbor- 
hood of the pair production region. Figure 6 shows that there exists a Vese- 
lago medium near the pair production region. Random points of normal and 
anomalous dispersion are found away from the pair production region. 

Figure 7 indicates the presence of a Veselago medium along with the 
normal dispersion of waves. Figure 8 represents most of the region admit- 
ting normal dispersion for the waves with high frequencies in a Veselago 
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medium. In Figure 9, anomalous dispersion of waves in a usual medium is 
found throughout the region whereas Figure 10 shows that the waves with 
higher angular frequency are normally dispersed in a usual medium. Figure 
11 represents a usual medium admitting anomalous dispersion of waves in 
most of the region. In Figure 12, small regions admitting normal dispersion 
of waves are found in usual medium. 

In Figure 13, a very small region shows normal dispersion of waves in 
usual medium. In Figure 14, presence of a Veselago medium is confirmed. 
Normal dispersion occurs for the waves with low angular frequency, anoma- 
lous otherwise. A Veselago medium admitting normal dispersion, except for 
low frequency waves, is found in Figure 15. In Figure 16, a Veselago medium 
exists throughout the region. Moreover, a small region near the pair produc- 
tion region admits anomalous dispersion. Most of the region shows normal 
dispersion. Figure 17 indicates anomalous dispersion of waves in a Veselago 
medium for most of the region. In Figure 18 random points of normal and 
anomalous dispersion are found in a Veselago medium. 

In Figure 19, a large and a small region exhibit the properties of a Vese- 
lago medium. Moreover, near the pair production region normal dispersion 
of waves exists. Figure 20 shows that most of the region near the pair produc- 
tion region admits normal dispersion of waves in a usual medium. In Figure 
21, random points of normal and anomalous dispersion are found. Further, 
small regions of a Veselago medium are present. In Figure 22, we obtain the 
usual medium near the pair production region. Random points of normal 
and anomalous dispersion are found. Figure 23 indicates that the medium 
is usual admitting anomalous dispersion of waves near the pair production 
region. Figure 24 shows anomalous dispersion of waves in the usual medium 
near the pair production region. 

From the Figure analysis, we arrive at the following conclusions. 

1. We can conclude that the magnetosphere allows the waves to disperse 
normally towards the event horizon as well as the outer end of the 
magnetosphere. Thus negative energy flux can fall inside the black 
hole and rotational energy can be extracted in response. 

2. It does not matter whether the region lies towards the event horizon or 
towards the outer end of the magnetosphere, the presence of a Veselago 
medium only depends on the rotation of the background. 

When we compare the results with the previous literature, we obtain the 
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following: 



• On comparing our results with [31j, only Figures 7 and 14 show that 
energy can be extracted by inflow of negative energy influx. In this 
work, we have obtained Figures 15, 7, 10 and 20 which show normal 
dispersion of waves. This indicates that there are more chances of waves 
to move towards the outer end of the magnetosphere (Figures 7, 10 and 
20). Consequently there are more chances of the extraction of energy. 

• Here we obtain clear cut conditions to show that extraction of energy is 
possible. In [32], only inflow of the waves towards the event horizon is 
possible and their outflow towards the outer end of the magnetosphere 
is not possible. Thus the exchange of energy from the event horizon 
with the accreting material was expected to be allowed. However, here 
we have obtained clear evidences that the energy can be extracted from 
the event horizon, passes through the accretion disk and can move 
towards the outer end of the magnetosphere as given by observational 
identifications in [23]- This does not negate the work given in [52] , 
but generalizes that idea. 
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Appendix A 



This Appendix includes the GRMHD equations and their perturbed form 
by using some assumptions. The component form of these equations is also 
given. 

The set of perfect GRMHD equations for isothermal plasma living in 
Kerr planar analogue (with rotation assumed to be along x-direction) can be 
written as (Eqs.(4.1)-(4.5) of [32]) 



^ + (B.V)/3 = Vx (VxB), 

dr 

V.B = 0, 

^^^ + (p+p)7^V.^ + (p + rtV.V = 0, 

Dt Ut 

iP + pW + ^] + (P + p)l'V,V^ - l^B^B, 



(Al) 
(A2) 

(A3) 



dV^ 
dr 



Hp + pWVi^kV^ + (p + ph^v.v^kV'' = (p + pW(^,,v^ 



-p, + -(i?,. 



An 



B X V X 



dB 

dr 



(A4) 



7'V.^(p + p) + 2{p + ph'Y.^ - ? + (P + P)7'(V.V) 
Ut Ut dT 

-(P + p)7'V.(V.V)/5 + ^ [(V X B).(V X B) 

+ (V X B).-^(V X B) + (V X B).{(V X B).V}/3 = 0, 

dr 



(A5) 



with V, B and E are fiducial observer (FIDO) measured fiuid velocity, mag- 
netic and electric fields respectively, ^ = ^ — /?.V is the FIDO measured 
rate of change of any three-dimensional vector in absolute space, 7 is the 
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Lorentz factor and = ^ + V.V = ^ + (V — /3).V is the time derivative 
moving along the fluid. 

The perturbed variables take the following form 



P = P + PP. 
V = V° + V, 



p = p +PP, 
B = B° + Bh, 



(A6) 



where the dimensionless perturbed quantities are 



P 

V 

b 



5p 



6p 



p{t,x,z), p= — =p{t,x,z), 
P P 
6Y = v^(t, X, z)e^ + v^{t, X, z)ez, 

— = z)ex + b^{t,x, z)e^. 



(A7) 



When we introduce the perturbations from Eq.( lA6p . the linearized GRMHD 
Eqs.(lAl])-(jA5l) become 

^ - /?.V)(5B)| = V X (v X B) + V X (V X 5B) - {SB.V)P, (A8) 

(A9) 



V.((5B) = 0, 
^ + (V - V 1^ {5p + 5p) + (p + p)7V. 



d_ 

di 



+ (V-/3).V V 



+(p + p)(V.v) + {5p + 5p)(V.V) + {5p + 5p)7V.(V.V)V 

= -2(p + p)72(V.v)(V.V) ln7 - (p + p)72(V.VV).v 

+ (p + p)v.VlnM, (AlO) 



(P + P)7' + 



B^ 

47r 



9_ 
at 



47r 



B X <^ V X 



- /?.VJ t;^ 

d{5B)\ 
dr J 



-(<^p),i + 7'[(<5p + <5p)V^' 



+2(p + j9)7'(V.v) V-'' + (p + , + —{B,, - B,,)5B^ 

4:71 



-ip + ph\v,V^+v,V')Vk,,V^-^^{i6p + 6p)V^ 

+2(p + p)^\V.y)V' + (p + pyjVi^j - i%{{6p + 6p)V' 

+4{p + p)J\Y.^r)V^ + {p + p^jV^^kV" , 



(All) 
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7' || + (V - /5).v| {5p + 5p) + 2(p + p)7^V. || 

+ (V-/5).V}v- (p + p)72v.Vlnu + (p + p)7^V.(v.V)V 

+6(p + j9)7^(V.v)V.(V.V)V + 2{5p + 5p)7^V.(V.V)V 

+2(p + j9)7S.(V.V)V - 2(p + p)7^(V.v)V.V Inw + 2(p + p)7^ 

x(V.v)(V.V) - + (5p + 5p)7'(V.V) + (p + p)72(v.v) 

-7'(/3.V)((5p + + 2(p + p)7^(V.v)(/?.VlnM) - 6(p + p) 
X76(V.v)V.(/3.V)V - 2(p + p)7V.(/5.V)V - 2((5p + (5p)7^V.(/3.V)V 
-(5p + 5p)72V.(V.V)/3 - (p + p)72V.(v.V)/? 

-2(p + p)7^(V.v)V.(V.V)/3 + [(v X B).(V X B) 

47r 

+ (V X (5B).(V X B) + (V X B).(V x 5B) 

0. (A12) 



+ (V X B). <^ — X B + V X ' 



c?r c/r J 

These equations show the change in the behavior of moving mass, velocity 
and magnetic field of the fiuid when the perturbations are involved. The 
component form of Eqs. flA8P - flA9P are given as follows. 

^ + Vb^^^ + ub^^^ = -u'b^ + i^'b^ + Va;,z - >^Vz,z - >^'vz, (A13) 
dr 

-r + ^^z,x + Ubz,z = \Vz,x - Vcc,x, (A 14) 

CLT 

bx,x + b,,, = 0, (A15) 

dp dp ^ ^ _ ^ ^ _ , ^ / ^ r / / 

p— + p— + pVp^^ + pVp^^ + pup^z + pup^z -[p-p){pu + pu 
dr dr 

+PM7^ + (p + ]9)7^ {^~^ ^ ^~Zt) ^^P^ P){7rfa;,a; + QVz,z 

+mV^7'K,. + Vz^,)} = -{p + p)^\[{'K + -f^V^)V' + 2^\Vu']v, 

v' 

+ (p + p)[(l - 27V)f?- - 27V\/V>„ (A16) 

dv dv 
{ (P + P)7'vr + ^}^ + {ip + ph^uV - ^A} ^ - 2ijuVbz,z 

+ { (p + ph^ir - ^} {Vv^,x + uv.^,,) + { (p + ph\V + i)\ ] 
X (yvz,x + uvz,z) + ^{(1 - V^)bz^^ - (1 - u^)bx,z} = -ipuu'bx 
+V'{A' + ud'}b^ - pp,^ - {pp + pp)7\{(l + 7V^)\/' + -i'^uVu'] 
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-(P + p)7^m{(1 + HV^)uu' + A-kVV'}v^ - [{p + p)7^ 
X [{(1 + 27V)(1 + 27V2) - 7V2}\/' 

+27^(1 + 27V)m1/m'] + ^mA>^, (A17) 
{ (P + P)7^f^ + ^A^} ^ + { (p + p)^\V - ^A} ^ 

+ { (P + P)7^^ - V^A^} (Vt;,,^ + Mt;^^^) + 2ip\uVh,^, 

+ {{p + p)i*uV + ^A} {Vv,,^, + ui;,,,) - ^A{(1 - V^)h,,, 

— (1 — 'U^)^^.^^} = — \uu')hx — ipXwd'bz — pp'j'^{uug 

+7 - Vf3'} - [pp^, + p'p + pp-f^{uu'g + 7 - V(3'}] 

-(P + Ph^[l^u\l + 4j^V^)V' - (1 + 27V)/3' 

+272My(l + 272u2)m>^ - [{p + p)72{-27\\//5' + ^(1 + A-f\^)u' 

+2-f\Q + -f^u'^)uVV'} - '^XuX']v„ (A18) 

7'p|^ + P(7' - + P7'{p'm + pu' + 2pM7V - P^^/5'} + 7'PP,X^ 

+'y'^upp^z + p^'^{up + m'p + 2pM7^(y9 — p7^mV/?'} + 'j'^pp^x'^ + 7'^upP,z 

+ ^{2(p + p)7V - ^x} + ^{2(p + P)l'u + ^Ax} 

+^x,.[(p + pWH + 27V^9} - + VxM'^ip + p)7V + ^x] 
+t;^,a;^?[2(p + p)7^M + V^Ax] + Vz,z{{p + p){q + l^u^) - i^Xux} 
+V'X{(1 - u^)b.,z - (1 - V^)K,x + 2uVbz,z)} 
+i}uK{\' - xu') + #.(-AV + ux^') + ^^.[(p + p)7'm{27V' 
+67V(^ - /3'(7r + 7V)} - ^A'] + t;.[^A'(A - ux) 
u' 

+ (p + pW{ + 27W + 67^mV + 7V - Vp'g}] = 0. (A19) 

u 

We have used the conservation law of rest-mass for three-dimensional hyper- 
surface given by Eg. (13.21) to simplify Eq. flA16l) . 



Appendix B 

Here we include the Mathematica program used to calculate a dispersion 
relation. This dispersion relation is obtained from the real part of the deter- 
minant. 

In the start, we must specify relative assumptions given in Section 3 to 
simplify the dispersion relations. In the following, we only consider the veloc- 
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ity components given by Eq. fl3.3p . We shall change the values of u, V, u', V 
to obtain the dispersion relations for the velocity components given by Eq. fl3.4p 
We have used the notations T, S, k, I for tanh(z), sec h{z), k^, kz 
respectively. The components gij are the coefficients of c/s in respective 
equations Eqs.(|i21), fl^ -M- 

p = l/2 
p' = D[p,z] 
p=l/2 
p' = D [p, z] 
I = -k 
T = m 

S = Sqrt[l — m^] 
/3=(T-1) 
/?' = 

u = l/3 (-T-(3 S^ + T^)) 

u' = 1/3 (-52 + (2 T)/(3 52 + T2)) 

V = l/3 (2 T- (3 S^ + T^)) 

V = 1/3 (2 + (2 r)/(3 + T^)) 
A = 1 

\' = D[\z] 

7= 1/(1-^2- 

Simplify [7] 



b = Det[ 



/ 911 


912 


913 


914 


915 


916 \ 


921 


922 


923 


92A 


925 


926 


931 


932 


933 


93A 


935 


936 


9ii 


fl'42 


9i3 


9u 


945 


946 


951 


952 


953 


954 


955 


956 


\ 961 


962 


963 


964 


965 


966 J 



Expand [%] 

The determinant of the coefficients of the constants ci, C2, C3, C4, C5 and cq 
from Eqs. (14.21) . (14. 3p . (I4.5p - (l4.8p gives a dispersion relation sextic in k^. This 
equation cannot be solved analytically for the exact solutions and is solved 
numerically. Here we have given commands used to solve the dispersion re- 
lation only for the region towards the event horizon. Similar commands with 
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the replacement of {z, —5, —1, 0.2} by {z, 1, 5, 0.2} can be used to calculate 
the roots for the region away from the event horizon. The commands are 
given as follows: 

R = NSolve[b== 0,k]; 

Table[{z, to, Part[R, 1]}, {z, -5, -1, 0.2}, {cu, 0, 10, 0.2}] 
Table[{z, lu, Part[R, 2]}, {z, -5, -1, 0.2}, {uj, 0, 10, 0.2}] 
Table[{z, u, Part[R, 3]}, {z, -5, -1, 0.2}, {cu, 0, 10, 0.2}] 
Table[{z, uj, Part[R, 4]}, {z, -5, -1, 0.2}, {a;, 0, 10, 0.2}] 
Table[{z, u, Part[R, 5]}, {z, -5, -1, 0.2}, {uj, 0, 10, 0.2}] 

Afterwards, each root is separated in the form of arrays for further ma- 
nipulation. Each root is approximated for each point of the two-dimensional 
meshes —5<z< —1, < a; < 10 and 1 < z < 5, < a; < 10, with equal 
step lengths 0.2 for z and uj. The numerical values in each region are used 
to approximate interpolation functions. 

The data for each root is converted and separately dealt in Mathematica 
by the following commands. The graph of the root admitting real values at 
each point of the mesh can be generated by the following procedure to obtain 
Figure 1 (e.g.). 

datapts = {{an, ai2, ...ai 51}, {^21, ^22, ■■■^2 51}, •••{a26 1, ^26 2, ■■■^26 51}}; 
r = Listlnterpolation [datapts, {{—5, —1}, {0, 10}}] 
p = Re[r[x, y]]; 
q = Im[r[x,y]]; 

pi = Plot3D[p,{z,-5,-l},{a;,0,10},AxesLabel^ {"z","a;"," "}] 
p2 = Plot3D[q, {z, -5, -1}, {uj,Q, 10}, AxesLabel ^ {"z","a;"," "}] 
p3 = Plot3D[^, {z, -5, -1}, {uj, 0, 10}, AxesLabel ^ {"z", "cu" , " " }] 
l = D[r[x,y],y]; 
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p4 = Plot3D[^, {z, -5, -1}, {cu, 0, 10}, AxesLabel ^ {"z" , "a;" , " "}] 



These commands generate x-components of the propagation and attenu- 
ation vectors, phase and group velocities for the specific set of data points 
in the form of graphs. These graphs are used to investigate the properties of 
the medium. 
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A B C D 




Figures 7-12 indicate the dispersion relations (related to the velocity comp- 
onents given by Eq.f l3.3l) ) in the neighborhood of the pair production region 
towards the outer end of the magnetosphere. 
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A B C D 




Figures 19-24 show the dispersion relations (related to the velocity compo- 
nents given by Eq. fl3.4p ) in the neighborhood of the pair production region 
towards the outer end of the magnetosphere. 
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